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1. Introduction 

Quantum field theories (QFT) in two space-time dimensions at finite temperature 
are an appealing and complex subject which has attracted much attention from both 
theoretical and experimental point of view. It could be useful to recall, that a finite 
temperature QFT can be obtained by considering the model on a cylinder with the 
radius of compactification proportional to the inverse of the physical temperature. 

Since the 2c? Ising model fuses together the simplicity and the full features of more 
complex models, we considered the most general, relevant perturbation of such model 
(both thermal and magnetic) in order to cope with the problem of finite temperature 
corrections to thermodynamic observables (in particular, as it will be clear later, such 
corrections are nothing but the finite size scaling behaviour of the model as a function 
of the radius of compactification). 

In this respect we recall that, taking advantage of integrability, the finite 
temperature study of the pure magnetic perturbation of the model was performed in [1] 
where an expression for the one point functions was proposed in the framework of Form 
Factors. Such results were found to be in agreement with the high precision numerical 
analysis of [2]. 

In this paper, we investigate the finite temperature features of the model in the 
region of the coupling constants (h representing the coupling to the spin operator, 
and t = \/3 — /3 c \//3 the deviation from the critical temperature) where the thermal 
perturbation is small with respect to magnetic one, t/h s ^ 15 — > 0. 

The main idea underlying our approach is to evaluate the Finite Size Scaling (FSS) 
on the cylindric geometry of the main thermodynamic observables (magnetization, free 
energy, magnetic susceptibility) by means of high precision numerical data coming from 
numerical diagonalization of the Transfer Matrix. Consequently, the finite temperature 
corrections to such observables, i.e. the FSS behaviour, are given by standard 
Thermodynamic Bethe Ansatz (TBA) prescriptions. 

Such an approach is known to work when the perturbations are separately switched 
on, giving rise to Integrable QFTs (see [1-3]). The contemporary presence of both of 
them destroys the integrability, however since we are in the region of the phase diagram 
where t/h 8 ^ 15 — > 0, we postulate that the functional forms coming from TBA still hold. 

We showed that this expectation is indeed correct (at least within the precision of 
our data), and as a side result we were able to evaluate the infinite volume quantities 
with higher precision than before [4]. It is important to point out that such functional 
form can be rigorously derived only when integrable theories are concerned, however 
we expect that when integrability is lost, they still hold under pair creation threshold 
which is just the regime explored in our simulations. 

The work is organized as follows: in sect. 2 we briefly present the model and the 
numerical transfer matrix technique; sect. 3 is devoted to the study of the FSS of some 
observables of the Ising model in magnetic field; in sect. 4 we present our numerical 
results; finally, in sect. 5 we give our conclusion. 
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2. The model and the technique 

The Ising model in a magnetic field at an arbitrary temperature is defined by the 
partition function 

Z((3,h e ) = e /3E (™.-> a " CTm+ ^ E " a ", (1) 

<7i = ±l 

where the spin variable a n takes the values ±1; the notation (n,m) represents nearest 
neighbor sites on the lattice; the sites are labeled by n = (n ,ni) and the two sizes 
of the square lattice are L and R (they are taken to be different because our transfer 
matrix calculations will treat the two directions asymmetrically); the total number of 
sites of the lattice will be denoted as iV = LR. 

The coupling /3 is the inverse of the temperature, while the magnetic perturbation is 
introduced by the coupling hi = Hf3, where H is the magnetic coupling. This model 
undergoes a second order phase transition when hi = and (3 reaches its critical value 

Pc [5] 

(3 C =\ \og(V2 + 1) = 0.4406868 .... (2) 

The observables which we shall consider are the following: 

the Magnetization 

M((3,h e ) = L^-QogZtfM)) = ^<5>>- ( 3 ) 



the Magnetic Susceptibility 

\ 

~ dh 



dM(P, h E ) 

x(P,hi) = wr • (4) 



the Free Energy 

f(P,h E ) = ^\ogZ(P,he)-f b . (5) 

where the bulk constant /& = 0.9296953982 . . . can be found in [5]. 

Since we will often refer to exact results coming from the QFT which describes the 
Ising model in the scaling limit, we recall that its action is given by 

A IM = A CFT + r J d 2 xe(x) + h J d 2 xa(x) (6) 

where the first term is the action of the Minimal Unitary CFT describing the critical 
point of the model [6] , h is the coupling constant associated to the magnetic field and r 
gives the deviation from the critical temperature (for a review on this subject see [7]). 
Such a theory is integrable either when h = 0, r ^ and one can prove that it describes 
a massive Majorana fermion, or in the case h ^ 0, r = which turns out to be a non- 
trivial theory of eight self-conjugated particles whose spectrum follows the numerology 
of the exceptional group E 8 as proved by Al. B. Zamolodchikov in [8]. Simple arguments 
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from CFTs are sufficient to conclude that in the general case /i ^ 0, r ^ integrabilityjj 
is lost [9]. 

In order to analyze the finite temperature behaviour of the model our choice was 
to employ the numerical transfer matrix technique. The basic idea is to rewrite the 
Boltzmann weight by means of the so-called transfer matrix T, which turns out to be a 
a positive symmetric 2 R x 2 R matrix such that 

Z((3,h e )= e /3E <".-) a " CTm+ ^ =trT L = X f ( 7 ) 

o"i=±l i 

where the Aj are the (real) eigenvalues of T (details can be found in [10-12]). The 
numerical computation of eigenvalues and eigenvectors of the transfer matrix enables 
us to compute all the observables we need, provided that we specify the values of hi, (3 
and R. As an example, the free energy 

f(P,h i ) = ^ I \ogZ((3,h i ) (8) 

can be written, in the limit L — > oo, as 

f((3,h e )^^\og\ max (9) 

which is exactly what we need for the numerical study of the thermodynamic observable 
as functions of the length of the cylinder R. 



3. Finite temperature results, the integrable case 

Taking advantage of the integrability of the model (jHJ) when (r = 0, h ^ 0), Delfino [1] 
applied Form Factor method to propose a compact expression for the finite size (i.e. 
finite temperature) corrections of the one point functions of the Ising model with 
pure magnetic perturbation (see also [13-15] for a discussion on the finite temperature 
analysis of 2d integrable Quantum Field Theory). 

Using these results, a given one point function ($)_r calculated on a cylinder of radius 
R is related to its infinite plane value ($)oo by 

ylr 1 = 1 + - E A* *o(r,) + 0(e- 2 -) (10) 

1 ' 1=1 

where r« = vrtiR and mi, with i = 1,2, 3, are the first three masses of Zamoldchikov's 
mass spectrum [8]. 

The main result of [1] is to compute exactly the universal constants Af . A numerical 
check of (jlOp has been performed by Caselle et al. [2] using Transfer Matrix technique to 
study the Ising model perturbed with a magnetic field. As already pointed out in [1], the 

t) We would like to stress that we will always refer to the integrability of the continuum theory because 
we will use exact results coming from the latter. We are not at all interested in the integrability of the 
lattice model since our approach is numeric and hence it is not related to the lattice integrability of 

CO- 
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one point function of the perturbing operator a can be obtained by TBA calculations. 
In fact, following TBA prediction for the FSS behavior of the free energy, it is possible 
to write down 

f(R) = /(oo) + ^^ + 0(e- 2 -) (11) 



, R 

i=i 



where E(ri) is given by 



E( n ) = i K x ( Ti ). (12) 

7T 

and Kx(r) is the modified Bessel function of the second kind. 
Hence, for the free energy we have 

which gives us the first few orders (under the lowest particle pair creation threshold) of 
the finite size corrections; the constants Ci are given by Zamolodchikov's mass spectrum 
(see ref. [8]). A first derivative with respect to the magnetic field gives us just expression 
fllOjl for the magnetization ff 

fir = 1 + E — W + A ° = & ( 14 ) 

Woo ri 7T 2Af 

then a further derivative gives us the functional form of the finite size corrections for 
the magnetic susceptibility 



X(OO) ri 7T 



/ 



The exact values of the constants are given in table ^ 



C x = 4.40490858 . . . 
C 2 = 7.12729179... 
C 3 = 8.76155605 . . . 
A f = 1.19773338... 



Table 1. Numerical values of the constants. 



The aim of the next section is twofold: on the one hand, we perform a numerical check 
for the finite size formulas for free energy and magnetic susceptibility in the case of the 
integrable model (for internal energy and magnetization this was proved in [2]); on the 
other hand we check that all the functional forms also hold in the case of the mixed 
perturbation considered in [4]. The latter case implies that the constants acquire a non- 
trivial dependence with respect to both he and t. Their functional form can be derived 
by the same scaling arguments employed in [4]. 

ff The constants Af coincide with the coefficients Af found by Delfino with the Form Factor approach 
in [1]. 
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4. Numerical analysis 

In order to check the behaviour of our data coming from the numerical diagonalization 
of the transfer matrix, we use the effective procedure presented in our previous work [4] 
and in [16]. 

First, for each value of he we fit the value of {&)r according to the functional form of 
the finite size corrections found in the previous section. As an example we can consider 
the susceptibility, hence the fitting function is given by (see previous section) 

*M = i _ V ^ [K ( n ) - 8 n K^n)] + 0(e~ 2ri ); (16) 

where we keep the quantities Af and x(°°) as f ree parameters. Our fitting procedure 
follows a very stringent list of acceptance criteria as presented in [4] . The first non trivial 
result is that the fits are acceptable (for all the observable considered here) even if the 
temperature is not set to its critical value. In this case the behaviour of the observables 
as function of R cannot be proved to hold rigorously because the corresponding QFT 
is not integrable. Notwithstanding this, in the regime t/h 8 ^ 15 — > considered here and 
within our numerical precision, such functional forms continue to hold. 

To the best of our knowledge the present is the first evidence about the validity of 
formulas like that of previous section for non-integrable theories. It could suggest that, 
when small perturbations of integrable QFTs are considered, then the finite temperature 
behaviour of the thermodynamic observables remains the same and the only price to 
pay is a change in the constants Af which acquire a non-trivial dependence with respect 
to the couplings t and he- The following quantitative results were found: 

(i) The numerical results for the infinite R observables /(oo), (er)oo, x(°°) are n °t 
only in complete agreement with respect to the estimates of [4], but they show an 
improvement in their precision of about one order of magnitude. As an example of 
this we quote the following results for the susceptibility 

= 8 r h = 0.07 



x(Pc, 0-07) = 0.7574405(3) (17) 
compared with the value obtained in [4] 

*(&,0.07) = 0.757441(1); (18) 
(3 = 0.4349 .. ., h = 0.07 

X(0.07, 0.4349 . . .) = 0.8465333(3) (19) 
compared with the value obtained in [4] 

X(0.07, 0.4349 . . .) = 0.846533(1). (20) 

We stress that it is remarkable that such an improvement comes independently of 
the fact that (3 is fixed at its critical value or not. Moreover such an improvement 
in precision could be expected since in this latter case we used the functional 
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form of the finite temperature corrections instead of fitting with a set of generic 
exponentially decaying functions. 

(ii) As anticipated, we supposed that the constants Af are functions of t and hi. In 
order to test such a conjecture, we first define the following subtracted quantities 
in order to discard all the corrections which depends only on the magnetic field hi, 

At(t)\ he - Af(0)\ he = J* (hi) t + J&ht) t 2 + J 3 %h e ) t 3 + . . . (21) 

As pointed out in [4] , it is a crucial step in order to compute the leading corrections 
due to the thermal perturbation. In this way we are able to fit our data with (|2*T|) . 
The precision of our data enable us to compute only the linear term of such an 
expansion. 

By means of the expansion for the scaling functions we are in the position to fit 
the amplitude J® {hi) as a function of the magnetic field alone. We recall that, in 
principle, the functional form of the <jfi(ht) is non-trivial to be obtained, however, 
they can be determined by the same Renormalization Group and CFT arguments 
given in [2, 4] in order to construct the scaling functions of the thermodynamic 
observables. The results are 

Free Energy: The fit with the scaling function 

Jiiih) = ^ + W 2 hf +W 3 hp + ... (22) 
hf 

and we can give a reliable estimate only of the leading term of the expansion 
Wi = 90(8). (23) 
Magnetization: The fit with the scaling function 

JlM = + K 2 hp + K 3 hf + . . . (24) 
hf 

gives us 

2.09 < H x < 2.15. (25) 
The other observables were outside our numerical precision. 

(iii) Finally, a further check about the validity is given by the analysis of the data at 
the critical temperature (3 C . In this case the scaling function for the coefficients Af 
depend only on the magnetic field hi and is given by 

Af(h e ) = Af(l + b lt *\h e \% + b 2 ^\h e \%) +.... (26) 

Employing them as fitting functions, we find 

8.07 < A\< 8.20 (27) 

16.12 < A{ < 16.28 (28) 

where the numerical values of the amplitudes perfectly agree with their theoretical 
counterparts: A\ = 8.09997. . . and A[ = 16.1999 . . . computed according to the 
formula? given in the previous section. Such results can be considered as further 
checks (beside the results of ref. [2] where magnetization and internal energy were 
considered) of the results coming from form factors. 
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Table 2. Value of the parameters used for the diagonalization of the Transfer Matrix. 



R 




X 


9 


0.; 


360965432368 


10 


0, 


352426578479 


11 


0, 


348936912437 


12 


0, 


347512060489 


13 


0, 


346931329088 


14 


0.; 


346695116084 


15 


0, 


346599363455 


16 


0, 


346560454648 


17 


0, 


346544860407 


18 


0.; 


346538199957 


19 


0.; 


346535213111 


20 


0.; 


346532714557 



Table 3. Magnetic susceptibility as function of the compactification radius L 
(7^ = 0.07,^ = 0.4349...). 

The complete list of the numerical values of the parameters we used for the 
numerical diagonalization of the Transfer Matrix is given in table El A significant 
sample of the data employed in the computation of finite temperature corrections is 
given in table El 

5. Conclusions 

Let us recall our main results. Our numerical study confirms the conjecture that 
the functional forms of the finite temperature corrections derived by means of 
Thermodynamic Bethe Ansatz for the Ising model with magnetic perturbation at t — 
(T = T c ) still holds when t ^ supposing that t/h 8 ^ 15 — > 0. We stress that, the 
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functional form of finite temperature corrections works well even if integrability is lost 
because we explored a region of the phase space where anelastic effect can be neglected. 
As explained in [17] the thermal perturbation results only in a correction to the mass 
spectrum of the integrable Ising field theory in magnetic field; no other effects due to 
the loss of integrability are expected. 

Finally, it could be worth to notice that, the high precision of our data allowed 
us to achieve the following numerical results: first we obtained the infinite volume 
estimates of the observables with typically one order of magnitude of precision more 
than the method previously employed in [4]; second we get information about the first 
subleading correction due to the reduced temperature. 
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